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FOURIER  INTEGRAL  METHODS  OF  PATTERN  ANALYSIS 


Abstract 

The  source  of  radiation  f(r)  and  Its  amplitude  pattern  g(l{)  are  a 
pair  of  Fourier  Transforms  j  f  (r)  g  »  where  g  (k)  =  J  f  (r)e****’‘dV 
and  f  (r )  OC  Jg  (k) e~^*^***do.  If  each  point  of  (r)  spreads  out  Into 
a  function  f«  (r) ,  the  resultant  source  Is  the  convolution  (Faltung), 
fi(r)<‘f,  (r),  while  Its  pattern  Is  the  product  gi(k)g2(k).  By  means 
of  this  theorem  the  solution  of  various  problems  can  be  written  down 
by  Inspection.  For  the  mattress  array,  ^0*^1  ***s*^»  fioSig»ga »  where 
Tq  Is  the  polarization  and  f 1 ,  f „ ,  fj  are  the  linear  arrays  along  the 
three  edges.  Other  problems  are:  (a)  Two  identical  antennas,  (b)  Bino¬ 
mial  distribution  of  points,  (c)  Uniform  distribution  of  points  (grat¬ 
ing)  ,  (d)  Trapezoidal  or  triangular  distribution,  (e)  Effect  of  a  linear 
phase  error.  Series  Expansions:  g(u)  =  2u,^u^/k!  where  u  =  (2n/X)sln  9 
and  Wic  •  x’*f(x)dx.  From  this  Is  derived  the  beam  width  W  (in 
radians)  for  n  db  down  on  the  main  lobe.  W  =  AVlid  -  nn)\/d  where 
A  »  0.152d  Vii„/  Mi ,  B  =  (1  -  p^,jj4/3ii2*) /34. 7  and  d  =  2a.  For  uniform 
lllumlnatlmn,  =  2/ (k+l)  and  W  *  0.629  Vnd  -  n/86.8)X/d,  correct 
to  ■  0.6%  In  W  down  to  26  db. 
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A.  INTRODUCTION 


1.  Purpose  and  Scope  of  Report 

The  development  of  high  gain  antennas  for  microwaves  has  brought 
Into  existence  problems  which  resemble  those  In  optics.  Although 
the  exact  solution  must  be  in  agreement  with  electromagnetic 
theory,  the  simpler  scalar  theory  of  physical  optics  or  sound  is 
sufficient  for  many  problems,  particularly  when  the  overall  dimen¬ 
sions  of  the  antenna  contain  several  wavelengths.  Since  the 
radiation  pattern  of  an  antenna  at  a  large  distance  can  be  ex¬ 
pressed  as  a  Fourier  Integral  of  the  currents  over  the  antenna.  It 
is  natural  that  certain  theorems  concerning  the  Fourier  Transform 
are  found  to  be  useful  In  yielding  short-cut  solutions.  These 
three  dimensional  theorems,  when  reduced  to  one  dimension,  are 
analogous  to  those  of  the  Operational  Calculus  which  has  long  been 
used  in  electrical  circuit  theory. 

The  present  papei^  Is  the  first  of  three  sections  on  Methods  of 
Diffraction  Pattern  Analysis  by  the  author  and  co-workers.  Follow¬ 
ing  an  elementary  Introduction,  several  Fourier  Integral  Theorems 
are  developed  which  form  the  basis  for  a  number  of  short-cut 
methods.  These  are  Illustrated  by  numerous  practical  examples 
which  include  the  effect  of  linear  phasing  and  polarization. 
Section  II  lists  tables  of  functions  which  are  particularly  useful 
in  calculating  the  diffraction  patterns  of  linear  arrays.*  Opera¬ 
tional  methods  are  extended  to  solve  for  the  effect  of  phase 
errors  across  the  aperture  on  the  diffraction  pattern.  Section  III 
deals  with  the  scattering  cross  section  of  various  shaped  objects. 

2.  Analogy  between  Microwave  Antennasand  Optical  Ref  lectors 

In  this  section  it  Is  assumed  that  the  reader  is  familiar  with 
geometric  optics  and  with  Huygens'  wave  theory  as  presented  In  the 
usual  college  textbook.  The  representation  of  phase  angle  cp  as  a 
rotating  vector,  e^'P,  is  a  useful  concept  In  studying  Interference 
and  diffraction,  while  the  concept  of  radiation  from  a  dipole  Is 
essential  In  any  electromagnetic  problem. 

Figures  1  to  4  illustrate  a  number  of  microwave  antennas  whose 
action  is  explained  by  elementary  optics.  In  Fig.  1,  rays  which 
start  from  the  focus  of  a  paraboloidal  reflector  are  reflected 
parallel  to  the  axis,  according  to  the  laws  of  geometric  optics. 
Thus,  the  use  of  paraboloidal  reflectors  In  flash  lights,  search 
lights,  reflecting  telescopes,  sound  detectors  and  reflectors  for 
microwaves.  Is  based  on  the  assumption  that  geometric  optics  Is 
valid.  Although  geometric  optics  Is  simple  to  apply.  Its  Justifi¬ 
cation  and  Its  limitations  depend  upon  the  more  fundamental 
Huygens'  wave  theory  which  is  illustrated  in  Fig.  2.  Huygens' 
Principle  also  explains  the  formation  of  a  cylindrical  wave  from 
an  array  of  point  sources,  (Fig,  3a)  and  the  formation  of  a  plane 
wave  from  a  broadside  array  (Fig,  3b) .  The  paraboloid  and  the 
parabolic  cylinder,  two  types  of  reflectors  commonly  used  in 
microwaves,  are  compared  in  Fig,  4,  In  Fig,  4b  a  linear  array  of 


TSi-l 


1 


Figure  1.  Paraboloidal  Fteflector. 
Law  of  Reflection 

Zr  "  Zi 


Figure  Paraboloidal  Reflector 
Optical  Path,  FA  +  iC  = 
DA  +  AT  =  constant 


Figure  3a.  Linear  Array  of  Point 
Sources  Form  Cylindrical 
Wave 


Figure"  3b.  Broadside  Array  of  Point 
Sources  Form  Plane  Wave 


Figure  4.  Comparison  of  (a)  Paraboloid  of  Revolution  and 
(b)  Parabolic  Cylinder. 
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dipoles  FF,  placed  at  the  focus  of  the  parabolic  cylinder,  emits 
a  cylindrical  wave  which  is  reflected  as  a  plane  wave. 

B.  POINT  SOURCES 


1.  Radiation  from  a  Point  Source 


The  intensity  of  a  wave  in  three  dimensional  space  at  a  large 
distance  R  from  a  small  source,  varies  according  to  the  usual 
inverse  square  law,  while  the  amplitude  varies  as  l/R.  An  actual 
source  may  be  a  collection  of  point  sources  or  an  equivalent  con¬ 
tinuous  distribution  of  point  sources  placed  at  a  reflector. 

In  the  case  of  longitudinal  waves  such  as  sound,  the  typical  point 
source  is  a  tiny  sphere  whose  radius  pulsates  so  as  to  emit  a 
spherical  wave.  Surfaces  of  constant  amplitude  <re  spheres  con¬ 
centric  with  the  source. 


In  the  case  of  electromagnetic  waves  which  are  transverse  waves, 
the  electric  field  E,  the  magnetic  field  H  and  the  direction  of 
propagation  are  perpendicular  to  each  other.  The  typical  point 
source  is  a  dipole  which  may  be  considered  to  be  a  very  short  wire 
of  length  L  carrying  a  sinusoidal  current  of  I  amperes.  At  dis¬ 
tances  large  compared  to  X/2ti,  where  X  is  the  wavelength,  the 
electric  field  in  volts  per  meter  is  given  by: 

E  =  60rt  IL  sin  w  (1) 

X  R 


\|/  is  the  angle  between  the  line  of  sight  and  the  dipole,  and  R 
is  the  range  in  meters.  The  ratio  L/X  is  assumed  small.  When  a 
collection  of  dipoles  have  parallel  axes  it  is  possible  to  treat 
them  first  as  point  sources.  Inserting  the  polarization  factor 
sin  at  the  end. 


The  motion  of  a  vibrating  sound  source,  or  the  current  in  a 
radiating  dipole,  may  be  represented  by  one  coordinate  of  a  point 
which  moves  In  a  circle  with  a  uniform  angular  velocity,  w  =  2Ttf, 
This  Is  represented  by  a  vector  which  has  turned  through  a  phase 
angle  qp  =  ut,  which  Is  equivalent  to  multiplication  by 


a.  Radiation  from  a  Collection  of  Point  Sources  ^ 


Fig,  B  Source  point  at  A,  Field  point  at  P, 
The  range  relative  to  the  origin  Is 
AR  =  r*fi. 
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In  Flg»  B  the  typical  point  source  A  is  a  distance  r  from  the 
origin  and  a  distance  R  from  the  field  point  P  which  in  Hirn  is  a 
distance  Rq  from  the  origin.  Since  the  radiation  travels  from  A 
to  P  in  the  time  R/c,  the  phase  at  P  Is  retarded  by  an  amount  (jR/c. 
The  field  at  P  due  to  a  number  of  such  point  sources  will  be,  ac¬ 
cording  to  Huygens'  Principle,  the  sum  of  the  Individual  fields^'® 


(2) 


If  we  let  R  = 
(2)  becomes 


The  external  factor  represents  a  spherical  wave  traveling  outward 
with  a  velocity  of  light,  since  the  phase  Is  constant  If  Rq  =  ct, 
and  the  amplitude  dlmlshes  as  l/Ro*  A  small  section  of  the 
spherical  wave  may  be  considered  part  of  a  plane  wave  with  a 
normal  n  In  the  direction  OP.  Henceforth  we  shall  omit  the  ex¬ 
ternal  factor  and  consider  only  the  contribution  within  the  sum¬ 
mation  sign,  which  Is  proportional  to  the  amplitude  of  the  re¬ 
sultant  plane  wave.  As  Ro  becomes  sufficiently  large  the  ratio 
AR/Ro  in  the  denominator  becomes  negligibly  small,  but  the  AR  In 
the  exponent  must  be  retained  as  the  phase  qp  relative  to  the 
origin  Is  not  small. 


Rq-AR  and  take  out  the  factor  common  to  all  terms. 


^  1-AR/Ro 


(3) 


qp  *  (i)  AR  =  8xfAR  “  271  ^  (4) 

c  c  X 

An  approximate  value  of  AR  Is  the  component  of  r  in  the  direction 
of  n  as  shown  In  Fig.  5 

AR  ■  r  cos  Y  ■  r-iT  (B) 

(r>n  Is  the  scalar  product  of  the  two  vectors  r  and  n,  while  y 
Is  the  angle  between  them) .  If  ^he  phase  change  per  unit  length  In 
the  direction  of  n  Is  designated  by  the  vector  propagation  con¬ 
stant  Hf  then  (4)  reduces  to 

qp  ■  2n:n*r  ■  lc*r  (6) 

and  the  field  (3)  Is  proportional  to  a  function  of  E. 

«<»«>  (7) 

Thus,  the  amplitude  of  the  equivalent  plane  wave  at  a  unit  dis¬ 
tance  from  the  origin  is  g(k),  which  is  the  sum  of  the  individual 
phase  vectors,  weighted  according  to  the  strengths  of  the  source 
points. 
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3-  Examples  of  Diffraction  Patterns ‘Derived  from  Point  Sources 

We  shall  now  consider  the  examples  of  Isotropic  point  sources 
Illustrated  in  Fig,  6.  By  definition,  a  unit  (isotropic)  point 
source  placed  at  the  origin  (Case  1),  radiates  a  unit  amplitude 
at  a  unit  distance  regardless  of  direction.  Hence,  g(k)  =  1. 
When  the  point  source  Is  displaced  a  distance  fi  from  the  origin 
as  in  Case  2,  the  phase  angle  Is  advanced  by  an  amount  q) 

=  (2nri /X)  COSY*  Note  that  <p  has  the  same  value  for  all  lines  of 
sight  that  make  an  angle  y  with  rj*  The  locus  of  these  directions 
forms  the  cone  Illustrated  In  Case  2. 

When  two  equal  point  sources  of  equal  sign  (Case  3)  ,  or  of  op¬ 
posite  sign  (Case  4) ,  are  placed  at  i  rj,  the  diffraction  patterns 
are  respectively  the  sum  and  difference  of  two  exponential  terms. 


Case  3.  g(V)»e  +e  =2cos  k»ri=2cos  q> 


(8) 


iTf-r,  -llf*?!  _ 

Case  4.  g  (k)  "e  -e  *21  sin  1cTi*21  sin  (p 


(9) 


One  must  remember  that  for  any  arrangement  of  point  sources  along 
a  line,  the  pattern  Is  identical  over  the  cone  of  directions 
shown  In  Case  2, 

The  diffraction  pattern  for  the  uniform  line  source  (Case  ») , 
extending  from  -ri  to  ri,  may  be  obtained  by  integrating  the 
patterns  of  a  uniform  density  of  point  sources  of  the  type  of 

Case  2. 


gOc) 


•/-ri 


lk*r 
e  dr 


1<P 

/.q, 


2  sin  (p 
|klco8  Y 


2ri 


cos  Y 
sin  (p 


(10) 


<P 


This  is  the  diffraction  pattern  of  a  single  silt  commonly  derived 
In  optics.  Other  examples  of  point  sources  will  be  taken  up  after 
the  discussion  of  the  properties  of  Fourier  Transforms. 

C,  THE  SOURCE  AND  ITS  PATTERN  AS  A  PAIR  OF  FOURIER  TRANSFORMS 
2.  General 

Wo  will  now  proceed  to  show  that  the  source  function  and  Its 
diffraction  pattern  are  Fourier  Transforms,  Let  us  designate  the 
source  function  by  f(r).  This  may  either  be  the  set  of  source 
points  described  above  or  a  continuous  function  f(r) 
of  the  position  vsctor  r.  The  amplitude  diffraction  pattern  (7) 
now  takes  the  form  of  a  Fourier  Integi'al 


‘The  term  diffraction  pattern  as  employed  in  optics  is  here  used 
interchangeably  with  the  term  radiation  pattern  which  is  employed 
in  radio. 
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<“(?)■  UnFfil 


Case  1.  Unit  Point  Source 
at  the  origin 


gW  -  ZAt, 


g(l^)  -  1 


Case  2«  Unit  Point  Source 
at  P  ■  ?! 


Case  3.  Two  Point  Sources, 
♦1  at  Pj ,  +1  at 


Case  4«  Two  Point  Sources, 
♦  1  at  Pj ,  -1  at-?! 


g(k)  »  =  e*<P 

(p  ■  kr^cosY  *  (arcTi/X) cosy 


g(k)  «  »  2  cos  9 


Zeros  for  q>  «  n/2,  3n/2,  etc. 


g(k)  ■  e'^‘^1  ■  21  sin  9 

Zeros  for  (p  ■  0*  n,  2n,  etc. 


Case  5.  Uniform  Line  Source 
from  -Pj  to  +?! 


Zeros  for  <p  =  n,  2Tt,  3ti,  etc. 


FIGURE  6.  DIFFRACTION  PATTERNS  OF  ISOTROPIC  POINT  SOURCES 
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g(lc)  =  y  f(r)e 


(11) 


IkTjiv 

The  Integration  Is  carried  out  over  the  volume  occupied  by  the 
source,  the  element  of  volume  being  dV, 

As  previously  mentioned,  g (k)  may  be  considered  to  be  the  ampli¬ 
tude  of  a  plane  wave  emitted  by  the  source  In  the  direction  of 
n  or  E.  If  we  reverse  these  waves  they  will,  by  the  principle  of 
reversibility  in  optics,  retrace  their  paths  exactly,  converging 
to  the  source  function  f (r) .  Thus,  we  may  represent  f(r)  as  the 
sum  of  plane  waves, 

f(r)  -  ^ g(k)e“^^*^do  ^2) 

The  integration  is  carried  out  over  the  surface  of  a  unit  sphere, 
the  element  of  solid  angle  being  do*  Note  that  the  k  In  the  phase 
factor  IS  reversed  because  of  the  reversal  of  direction.  We  have 
here  ignornd  a  normalization  factor  on  the  right  side  of  (12) • 

Equations  (11)  and  (12)  are  symmetric  In  their  mathematical  form 
and  are  known  as  a  pair  of  Fourier  Trans -forms.  This  fact  Is  some¬ 
times  denoted  by  the  double  arrow  symbol 

t(r)^  g(ia  (13) 

Again,  the  statement  that  "  g(E)  Is  the  Fourier  Transform  of  f (?) ” 
may  be  denoted  by 

g(k)  *  T  f(r)  (14) 

a.  Fourier  Transforms  in  Furious  Coordinote  Systems 

When  the  amplitude  function  is  known  over  a  simple  geometric 
surface  such  as  a  sphere  or  a  cylinder.  It  can  be  expanded  In  a 
series  of  orthogonal  functions,  each  of  which  has  Its  own  Fourier 
Transform.  These  are  given  by  Stratton*. 

The  case  of  rectangular  coordinates  is  quite  simple.  Let  the 
amplitudes  of  the  1,  J  and  k  components  of  the  vector  E  be  u, 

V  and  w.  If  x,  y  and  z  are  the  corresponding  components  of  r, 
the  phase  becomes 

V  ■  E*r  “  (ul+vj+wk)  •  (xl+yj+zk)  =  (ux+vy+wz)  (16) 

The  volume  element  dV  is  dxdydz  while  the  volume  element  in 
diffraction  space,  assuming  a  variation  In  wave  length,  becomes 
dudvdw/ (2ti) * •  If  we  divide  the  numerical  factor  (2n)~"  Into  two 
symmetrical  factors  (2tc)**"'^  we  arrive  at  the  usual  three  dimen¬ 
sional  transforms* 

g(u,v,w)  -  ^ fff  f(x,y,z)e^^’“*'^y**'*^ dxdydz  (16) 
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f(x,y,z)  = 


1 


dudvdw 


(17) 


[2n)  9/* 


/// 


In  case  f(x,y,z)  has  zero  thickness  In  the  i,J  ur  k  direction 
g(\i,v,w;  will  be  independent  of  the  corresponding  u,v  or  w  vari¬ 
able  due  to  the  fact  that  this  variable  vanishes  from  the  ex¬ 
pression  for  the  phase,  Uf>) .  Thus,  in  the  case  of  two  dimensions, 
which  Is  particularly  useful  for  broadside  arrays 


g (u,v) 

=  —  ^  f(x,y)ei<’«*vy>dxdy 

2n 

(IS) 

f  (x  ,y) 

=  —  yy  g  (n,v)  e~^  dudv . 

(19) 

In  the  case  of  one  dimension,  used  in  the 
arrays 


f  (x)  *  ^==  y^g(u)e~^’“du  . 


f  f(x)e^“*dx 


analysis  of  linear 

(20) 

(21) 


Usually  in  antenna  problems  we  wish  to  keep  the  frequency  constant 
which  Imposes  the  condition  that 

S  2  2  2 

u  +  V  +  w  =  k.  (22) 

This  condition  restricts  the  diffraction  space  to  the  surface  of 
the  sphere  mentioned  in  connection  with  the  Inverse  transform  (12). 

Again,  we  may  wish  to  keep  the  direction  fixed  and  merely  vary  the 
frequency*  In  case  the  direction  Is  along  the  x-axls,  the  one 
dimensional  problem  reduces  to  the  frequency  spectrum ^ g (u)  of 
a  transient  f(t)  as  treated  In  electrical  circuit  theory  •  There 
is  a  change  of  scale  arising  from  the  relationships  (p  ”  ux  ■  ut, 
X  ■  ct. 


g(u)  ^  yf(t)el“^dt  (23) 

f(t)  ,  j==:  /g(«)e“*^'^^du  (24) 

t#  Sft  ^ 


The  normalization  constants  for  any  of  the  above  equations  are 
chosen  so  that  there  Is  conservation  of  energy  or  conservation  of 
power  in  both  the  source  space  and  diffraction  space. 
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3.  Effect  of  Scaling 

If  the  amplitude  of  a  function  f(x)  Is  multiplied  by  the  factor 
8  while  it  is  made  narrower  in  the  same  ratio  so  as  to  preserve 
its  area,  its  Fourier  Transform  will  be  simply  broadened  by  the 
factor  Sf  while  its  peak  is  unchanged.  In  mathematical  notation, 
let 


sf(sx) 


g(u) 

‘  g(u/s) 


To  prove  this,  let  sx  =  x*  and  u/s  =  u’,  then  ux  “  u'x'.  The 
Fourier  Transform  of  sf(8x)  is 

_i _  f  8f(8x)e^^dx  =>  ^.1 _  f  f (x')e^“**'dx'  ■  g(u*) 

iTsr 

which  proves  (27).  This  principle  may  be  applied  to  any  of  the 
three  dimensions  independently. 

4.  Fourier  Transforms  of  Even  and  Odd  Functions 

Since  cos  ux  ♦  1  sin  ux,  g(u)  can  bo  written  In  the  form 


g  (u) 


fsirr- 


(x)  cos  ux  dx 


+  i  f(x)  sin  ux  dxj 


ge(u)  ♦  igg  (u). 


The  first  integral  Is  cero  for  an  odd  function  and  therefore 
applies  only  to  an  even  function,  f^,  while  the  second  applies 
only  to  an  odd  function,  f..  The  presence  of  the  1  indicates  that 
the  two  transforms  are  90^  out  of  phase.  This  was  discussed  in 
■ore  detail^"  In  R.L.  Report  T-7.  It  follows  that  the  power 
pattern  is  the  sub  of  squares  of  the  individual  patterns. 

P(«)  ■  g(u)  g(u)  -  (ge+lg„)  <go-Uo>  “  «e**  80* 

This  is  a  real  positive  function  with  no  zeros  since  g^  and  g^ 
are  unlikely  to  have  coincident  zeros. 

It  will  be  the  policy  In  the  remainder  of  this  paper  to  use 
particular  coordinate  systeas  only  when  they  simplify  the  problem. 
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5-  Fourier  Transform  of  the  Complex  Conjugate 

The  power  P(k)  In  the  diffraction  pattern  is  of  practical  im¬ 
portance  in  antenna  theory  and  measurements.  This  is  the  square  of 
the  absolute  magnitude  of  g (1)  or  the  product  of  its  complex 
con  jugates . 

pOf)  >  |g(t)l’  =  ga)g(i)  (30) 

Now,  to  find  the  complex  conjugate  g  (Tc)  we  simply  replace  the 
quantities  in  (11)  by  their  complex  conjugates  which  Includes 
reversing  the  i  in  the  exponent.  Thus, 

g(Tf)  =/f  (r)e"^^  (31) 

We  may  make  the  exponent  positive  again  by  reversing  all  k's  or 
reversing  all  r's.  These  changes  lead  respectively  to  the  trans¬ 
forms 

f(?)  g(-l) 

?(-r)i=t;f(Tc)  (32) 

6.  Effect  of  a  Linear  Phase  Change 

Let  g  (It)  be  the  diffraction  pattern  when  the  source  function 
f 1  (?)  is  real,  i.e.  all  source  points  are  vibrating  in  phase.  The 
effect,  then,  of  exciting  them  by  means  of  a  plane  wave  of  wave¬ 
length  Xg  traveling  in  the  direction  n^,  is  to  Introduce  a  linear 
phase  lag  (p^  at  each  point  (See  Fig.  7). 


R 


Fig,  7.  Path  Difference  when  Excited  by  Plane  Wave 
OC-OB  =  r.n-r»nQ=  r«(n-nQ)  =  r»N 


OB  2tc  _  _ 

qPo  =  2n-^  =  ~  r-n^j  =  r^k^, 

We  then  have  a  complex  source  function 
f(r)  =  e"^*”*‘ofi  (?) 


(33) 


(34) 
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j.eh  results  in  the  diffraction  pattern 

g  (0  =/f  (r)  e^'*  ®  dV  .yi ,  (?)  dV  -  (2_S^)  (M) 


By  adding  Eg  to  each  abscissa,  ve  obtain 


g(t*Eo>  “ 


(36) 


Thus,  the  effect  of  a  linear  phase  error  Is  to  shift  each  point  In 
^e  diffraction  pattern  gjd)  fron  the  direction  associated  slth 
Jc  to  that  of  (fc+kg)  •  This  theorem  for  three  dimensional  space 
corresponds  to  the  Heaviside  shift  theorem  in  the  operational 
calculus  of  electric  circuits. 


In  the  case  of  a  linear  array  excited  by  a  wave  traveling  along 
the  wave  guide,  the  phase  9  Is  shifted  by  an  amount  Brcr^/Xg, 
where  X.  is  the  wave  guide  wavelength.  Thus,  the  resultant 
diffraction  pattern  Is  a  new  function  of  ip. 


g(<p)  -  gi(<P-(Pg) 


(37) 


where  g, (9)  is  the  pattern  for  an  array  In  phase.  (We  are  here 
using  g(q^)  and  g(k)  Interchangeably  although  they  really  are 
different  due  to  the  scale  factor  between  <p  and  E) .  Since  g^  ((p) 
has  a  maximum  when  qi  Is  zero,  g((p-<pg)  will  have  a  maximum  when 
(<P-Tg)  “  0,  or  when 


cos  Y 


(38) 


V  Is  the  velocity  In  the  guide.  This  Is  Illustrated  In  Fig.  8a 
where  a  wavelet  travels  a  distance  OA  In  the  guide  while  It  is 
traveling  a  distcuice  OB  outside.  The  line  AB  Is  a  common  wavefront 
for  all  wavelets  along  the  guide.  The  same  constructlwi  Is  used  to 
explain  the  critical  angle  at  the  boundary  of  a  dielectric, 

bow  waves  from  boats  on  the  water,  or  the  shock  waves  from  bullets 
the  air.  Fig.  8c. 


(b)  CRITICAL  ANGLE  IN  GLASS 


(c)B0W  WAVES  FROM 
BOAT  OR  BULLET 


^ig.  8.  Huygens'  Have  Construction  for  Linear  Phase  Brrors 
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If  X  ==  Xgf  then  the  peak  intensity  occurs  when  y  Is  zero  and  the 
result  is  an  end-fire  array.  An  interesting  example  of  an  end  fire 
Is  the  case  of  two  point  sources  with  a  spacing,  d  >  X/4  and  with 
a  phase  difference  of  n/2,  or  1  n/4  from  that  of  the  origin.  The 
solution  Is  to  shift  the  curve  gj  (9)  >  cos  9  to  the  right  an 
amount  9.  *  k/4,  so  that  the  peak  value  occurs  at  y  *  0,  and  the 
first  and  only  zero  at  y  *  180®. 

D.  THB  CONVOLUTION  TRBORBM 

1.  Derivation  of  fi*ff5:±  gaSt 

The  convolution  theorem,  known  in  mathematics  as  the  Parseval 
theorem.  Is  a  powerful  tool  for  solving  many  diffraction  problems. 
The  following  ^derivation  for  three  dimensional  space  holds  equally 
well  for  one  or  two  dimensions.  Let  us  consider  the  two  source 
functions  ft  (r)  and  ft  (i)  of  Fig.  9. 

fi(r)  -  [a„r„]  ,  f ,  (i)  -  [b„i„]  (39) 

Thus,  f 1  (r)  Is  a  "set*  of  source  points  with  weights  and 
position  vectors  r^  In  three  dimensional  space,  while  ft (i)  has 
wslghts  b„  and  position  vectors  a^.  Note  that  the  square  bracket 
Is  employed  for  a  "set*  instead  of  the  1  which  would  have  denoted 
a  sum. 

The  convolution  of  fa  and  ft,  also  termed  the  resultant  or  the 
"Faltung*  (folding),  consists  of  spreading  each  point  of  one 
of  the  functions  out  into  the  other  function.  This  is  illustrated 
in  Fig.  10a.  Thus,  the  point  a„r„  is  spread  into  the  series 
•nC^m<V*m>l  ^^ile  the  convolution  f  is  given  by 

f  ■  f.*f.  •  K  («' 

The  star  between  two  fVinctions  indicates  their  convolution.  This 
is  symmetrical  with  respect  to  fa  and  ft,  the  typical  term  having 
a  weight  a|^b|g  and  a  position  vector  (^n'^^m^  '  Fig.  10b). 
Substitution  of  these  values  in  (7)  yields  the  diffraction  pattern 


g 


*  ^  *n  H  ^’m* 


ik.(r„+^) 


(41) 


But  this  is  identical  with  the  product  of  the  transforms  of 
ft  and  ft. 


8t«t  "  (43) 

.  i1c«(rn'''8_\  _ 

Since  in  either  case  the  typical  term  is  •  Thus 

we  arrive  at  the  convolution  theorem  which  states  that  the  trans¬ 
form  of  the  convolution  is  the  product  of  the  transforms . 


ft 


«a«. 


(43) 


12 


78B-1 


Figure  9. 


Point  functions  (?)  and  (s) ,  where  r 
position  vectors  in  three  dimensional 


and  s  are 
space . 


Figure  loa.  Convolution  of  f^  and  f^ 

f  =  f,*f. 


Figure  10b.  Typical  term  In  convolution. 

Weighted  points  a„b^  at 
position  +  i^. 


FIGURES  9,  10.  CONVOLUTION  OF  POINT  SOURCES  IN  SPACE 
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This  theorem  holds  equally  well  for  continuous  functions*  It  Is 
quite  apparent  that  we  may  form  the  convolution  of  a  third  func¬ 
tion  fa  With  f,  •'/hich  is  already  the  convolution  of  fi  and  fj* 
Again,  the  order  of  operations  is  Immaterial. 

*’i’  fa  .  .  .  •*-  •  ^  gigaga*  •  ‘gn 

Due  to  the  symmetry  between  the  f  and  g  functions,  the  above 
convolution  properties  hold  If  we  interchange  the  f's  and  g's. 

fift  52:  gi*ga  (46) 

The  convolution  theorem  Is  so  simple  that  frequently  by  the  time 
one  has  clearly  stated  the  problem  the  solution  can  be  written 
down.  Numerous  illustrations  of  the  method  are  given  later  In 
the  paper. 

a.  Integral  Representation  of  the  Convolution 

We  shall  now  derive  a  mathematical  expression  for  the  convolution. 
Let  us  change  notation,  replacing  r  by  r, ,  i  by  Fa ,  letting  r 
stand  for  the  vector  sum  (ri-tFa).  We  now  allow  the  volume 
element  fj (ri)dVa,  considered  as  a  mass  point,  to  spread  out  into 
the  function  fa(rg)*  The  convolution  f (r)  is  the  integral  over 
all  such  elements  of  volume* 

f(r)  =  fi(ri)  *  taGa)  =■  ff  i  (rt)  f  s  (^a )  dV^ 

(46) 

(fi)f*  (r-fj)  dVi 


If  the  functions  were  constrained  to  one  dimension  we  would  have 
the  usual  expression  for  the  convolution  of  f ,  (x)  and  fa (x) 

f(x)  =  f  1  (x)  *fa  (x)  =f  f  1  (a)  fa  (x-a)da  (47) 

where  a  is  a  variable  of  integration. 

3.  The  Parseval  Formulas.  Conservation  of  Energy 
Writing  (46)  as  the  transform  of  giga 

^fi  (F,)  fa  (F-Fi)dVi  =y^gi  (E)ga  (ic)e“^^*^dO  (48) 

If  )rj  =  0 

y*f  1  (Fi)  fa  (-Fi)  dVi  = /gi  (E)  ga  (K)  do 

Now  let  fa  (-Fi)  =  ^i(-r,),  whereupon  ga  (ic)  becomes  gi  (E)  according 
to  (32) . 
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Dropping  the  subscript 

/f(r)f(r)dp  =  /g(Ic)g(K)(l<r 
or 

J|f(T)l“  dp  =  /)g(k)|“  do-  (49) 

The  above  Parseval  formulas  are  usually  developed  in  the  lit¬ 
erature  for  one  dimension  only.  Eq,  (49)  is  really  a  statement  of 
the  conservation  of  energy  between  the  source  space  and  diffrac¬ 
tion  space,  or  between  antenna  primary  patterns  and  secondary 
patterns. 

4.  The  Power  Diffraction  Pattern  as  the  Transform  of  f(r)*?(-r) 

In  most  measurements  and  experiments  involving  diffraction  pat¬ 
terns,  the  phase  of  the  pattern  is  not  detected,  only  the  power 
P(k)  =  g(i:)g(k)  being  significant.  But  according  to  (30),  (32)^  and 
(43),  PCS)  is  the  transform  of  the  convolution  of  f(r)  and  f(-r) 
showing  that  this  is  the  significant  characteristic  of  the  antenna 

f  (?)*?(-?)«  g(E)g(Ic)  =  P(l£)  (00) 

The  same  equation  is  obtained  from  the  Parseval  Formulas  (46,  48) 
by  stibstituting  g(Ic)  for  gi  (k)  and  g(k)  for  ga  (t)  • 

E.  EXAMPLES  OF  CONVOLUTIONS  IN  SPACE 
i.  Effect  of  a  Displacement 

The  effect  of  a  displacement  of  the  origin  of  a  source  function 
fi (r)  to  a  position  r^  can  be  considered  to  be  the  convolution  of 
and  f  1  (r) 

f(r)  =  f a  (r-r^j)  =  r^  f  1  (?)  (61) 

The  resultant  diffraction  pattern  is  then  the  product  of  the  two 
patterns 

g(E)  =  g^(5)  (62) 

Note  that  this  Introduces  a  Linear  phase  error  into  the  diffrac¬ 
tion  pattern. 

s.  Effect  of  a  Linear  Phase  Error 

This  effect,  discussed  earlier,  is  the  same  as  the  above  example 
with  the  roles  of  f  and  g,  and  those  of  k  and  r  reversed.  The 
resultant  source  function  Is  the  product 

f(r)  “  fj(f) 

while  the  resultant  diffraction  pattern  is  that  derived  earlier 
by  a  direct  method.  (See  (36). 

g(Ic)  =  gi  (K-Eq)  (e2a) 
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3.  Effect  of  Polarization 

If  the  currents  in  an  antenna  are  all  piirallel  we  may  consider 
that  each  scalar  point  has  been  spread  out  Into  a  dipole*  The 
resultant  source  function  is  the  convolution  of  the  scalar  source 
functlofk  f(r)  and  the  dipole,  so  that  the  diffraction  pattern  Is 
the  product  of  the  scalar  diffraction  pattern  and  the  factor 
sin  \|r  of  (1) . 

4,  Broadside  Mattress  Type  Array 

The  mattress  type  array  of  Fig.  3b  consists  of  two  rectangular 
layers  of  equally  spaced  dipoles  with  their  axes  parallel  and  with 
corresponding  dipoles  in  the  front  and  back  layers  arranged  for 
end-fire.  If  all  dipoles  are  equally  excited,  the  mattress  Is  the 
convolution  of  four  functions,  the  polarisation  f^  and  the  three 
linear  arrays  along  the  three  edges.  The  diffraction  pattern  Is 
therefore  the  product  of  the  four  patterns  according  to  (44). 

1*0  fa  fa  SIS  (53) 

This  equation  applies  to  a  much  wider  class  of  antennas  than  the 
conventional  mattress  type,  since  the  individual  axes  may  have  any 
orientation,  and  the  spacing,  amplitude  and  phase  along  any  one 
linear  array  may  be  arbitrary*  The  ease  of  a  three  dimensional 
grating  with  inclined  axes  is  treated  by  Stratton* 

5.  Tuo  Identical  Antennas 

The  arrangement  of  two  identical  antennas  placed  with  centers  at 
±  ra  is  simply  the  convolution  of  one  of  the  antennas  with  two 
points  located  at  ±  Tg *  The  diffraction  pattern  is  then  the 
product  of  the  pattern  gi  (Tc)  for  one  of  the  antennas  and  the 
Interference  pattern  for  two  point  sources,  either  In  phase,  as 
in  Case  4.,  (Fig.  6), or  with  an  arbitrary  phase  difference  due 
to  the  length  of  feed  line.  If  we  also  include  the  polarization 
factor  we  can  write  the  complete  pattern  as  the  product  of  three 
factors 

g(0  “  gi  (®  cos(lc-ra-<Po)  sin  (54) 

The  next  paragraph  discusses  a  special  case. 

6,  Double  Slit  Source  Considered  as  a  Convolution 

The  double  slit  source  of  optics,  illustrated  in  Fig.  11,  can  be 
considered  to  be  the  convolution  of  a  pair  of  points  situated  at 
±a  and  a  single  slit  of  width  2b.  The  over-all  diffraction  pattern 
Is  the  product  of  the  two  patterns  2  cos  au  and  2  (sin  bu) /u 
illustrated  in  Fig.  6.  It  Is  recalled  that  the  first  factor  is 
unity  when  au  is  a  multiple  of  n  or  when  the  path  difference 
(2a  sin  0)  ’°(2a  cos  centers  Is  a  whole  number  of  wave¬ 

lengths,  so  that  the  two  sources  are  in  phase.  The  second  is  zero 
when  bu  is  e  multiple  of  n  or  when  the  variation  in  path  length 
(2b  sin  0)  over  a  single  slit  is  a  whole  number  of  wavelengths. 
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Diffraction  Patterns  For: 

Two  Points,  gi  (u)=2  cos  au 
One  Slit,  (u) =2  sin  bu 

u 

Two  Slits  g(u)=  gi  ga 

“  4  cos  au  Sin  bu 
u 


Fig.  11.  Double  Slit  as  the  Convolution 
of  Two  Points  and  a  Slit 


7.  Double  Slit  as  the  Difference  of  Two  Slits 

An  alternate  method  of  attack  is  to  make  use  of  the  fact  that 
Fourier  Transforms  are  linear  transformations,  i*e.  the  transform 
of  a  sum  of  functions  is  the  sum  of  their  transforms. 

fi  (x)  +  fa  (x)  +..«  giOO  +•••  g2(u)  +•••  (86) 

Fig.  12  shows  how  the  difference  of  two  slits  of  widths  2(a'*’b)  and 
2(a-b)  is  equivalent  to  the  double  slit  of  Fig.  11.  The  corres¬ 
ponding  diffraction  patterns  in  the  two  figures  are  also  equiv¬ 
alent  since  they  are  trigonometric  identities. 


g(u)  _  2  sin  (a^b)  u  2  sin(a-b)u 

u  u 

4  cos  au  sinbu 

U 


Fig.  12.  Double  Slit  as  the  Difference  of  Two 
Slits  of  Widths  2(a'*’b)  and  2(a-b) 
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F.  CONVOLUTIONS  OF  PAIRS  OF  POINTS  IN  ONE  DIMENSION 

1,  Successive  Convolutions  of  Pairs  of  Points 

Fig.  13  illustrates  successive  convolutions  of  pairs  of  points  In 
one  dimension.  The  top  row  shows  one  pair  at  ±a  with  a  diffraction 
pattern  cos  au.  The  next  row  shows  Its  convolution  with  a  second 
pair  of  points  at  ±  F,  the  resulting  pattern^belng  cos  au  cos  bu. 
The  third  row  Is  the  convolution  of  the  second  row  with  a  pair  of 
points  at  1  Cf  the  pattern  being  cos  au  cos  bu  cos  cu.  This 
process  may  be  continued  Indefinitely. 

The  final  pattern  Is  zero  whenever  one  of  the  component  patterns 
Is  zero.  If  a,  F  and  c  differ  slightly,  their  first  zeros  differ 
slightly  so  that  the  pattern  Is  almost  zero  over  a  range  of 
angles.  The  pattern  can  become  unity  at  a  point  other  than  the 
origin  only  If  a,  S  and  c  have  simple  Integral  ratios.  Thus,  the 
presence  of  side  lobes  due  to  equidistant  supports  can  be  reduced 
only  by  a  careful  study  of  vaxlable  spacing. 

3.  Binomial  Distribution  of  Points 

Fig.  14  shows  the  binomial  distribution  of  n  points  resulting 
from  (n~l)  successive  convolutions  of  a  pair  of  points  with 
Itself.  Thus,  a  b  =  c  etc.  The  points  are  equally  spaced  and 
their  amplitudes  are  the  binomial  coefficients  for  terms  In  the 
expansion  of  (l+l)”/2''.  The  diffraction  pattern  Is  simply  cos”au. 
The  position  of  the  zeros  and  of  the  points  where  the  pattern  is 
unity  Is  the  same  regardless  of  n.  If  the  spacing  2a  equals  or 
exceeds  X,  strong  higher  orders  are  formed  as  in  the  case  of  the 
diffraction  grating.  When  2a  Is  less  than  X/2  the  diffraction 
pattern  Is  a  monotonlc  decreasing  function  approximating  a  Gauss 
Error  curve,  but  the  array  Is  too  weakly  illuminated  near  the 
ends  to  be  practical. 

3.  Grating  of  2*'  Points  (Fig.  15) 

Again,  if  b  =■  2a,  c  =  2b,  etc.  we  find  that  (n-1)  convolutions 
produce  a  diffraction  grating  of  2”  points  of  equal  amplitude 
equally  spaced.  The  resulting  pattenn  is 

g(u)  °  cos  au  cos  2au  cos  4au  etc.  (56) 

Since  a  high  density  of  points  Is  equivalent  to  a  uniform  line 
source,  It  follows  that  the  pattern  for  a  uniform  line  source  of 
width  2w  can  be  expressed  as  an  Infinite  continued  product. 

sin  wu  cos  wu  cos  wu  cos  wu 

•  •  • 

u  2  4b 

4.  Grating  of  n  Equal  Points  as  a  Convolution  (Fig.  16) 

1  3 

A  simple  derivation  of  the  diffraction  grating  formula  is  obtain¬ 
ed  as  follows.  Let  ^^(x)  designate  an  array  of  n  equidistant  point 
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Types  of  f(x) 


Figure  13  •  Successive  ^convolutions  of  pairs  of  Points 
Separations  ya,  2b,  2c,  etc. 


cos  nu 


7 

c  os  au 


3 

cos  au 


Figure  14*  Binomial  Distribution 

Separations  a  =  b  =  c,  etc. 


Figure  15-  Uniform  Grating  of  2”  Points 

Separations  b  =  2a,  c  =  2b,  etc. 
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au 
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Figure  10-  Diffraction  Grating  of  n  Points 


FIGURES  13-16.  CONVOLUTIONS  OF  PAIRS  OF  POINTS 


sources  of  constant  amplitude  separated  by  Intervals  2a.  Then  let 
each  point  be  spread  Into  a  slit  of  width  2a.  The  result  is  one 
long  slit  of  width  2w  =  2na>  the  diffraction  pattern  of  which  is 
(sin  wu)/u  =  (sin  nau)/u.  But  the  long  slit  was  formed  from  the 
convolution  of  the  point  pattern  and  the  single  slit  of 

width  2a,  so  that  the  resultant  pattern  can  also  be  expressed  as 
the  product  of  the  point  pattern  g^(u)  and  the  single  slit  pattern 
(sin  au)/u.  On  equating  the  two  expressions, 

g,,(u)  sin  au  _  sin  nau 

u  ■  u  (87) 

Hence 

gjj(u)  _  sin  wu  _  sin  nau  (88) 

sin  au  sin  au 

This  reduces  to  n  when  au  is  zero  or  a  multiple  of  n,  l.e.  when 
waves  from  adjacent  openings  differ  by  a  whole  number  of  waves. 
The  pattern  is  zero  whenever  wu  is  a  multiple  of  n.  This  causes 
(n-2)  side  lobes  to  appear  between  each  of  the  strong  peaks. 
(See  any  optics  text  for  illustrations).  Midway  between  the  strong 
peaks,  the  envelope  of  g,j(u),  which  is  l/(sln  au)  ,  is  unity  so 
that  the  minimum  side  lobe  has  a  relative  amplitude  of  1/n  or  a 
relative  intensity  of  (l/n)  •  The  corresponding  intensity  pattern 
from  the  slit  of  width  2w  would  have  been  (2/n) ’  as  much,  or 
4  db  lower. 

5*  Grating  of  n  Equal  Points  as  a  Product 

The  diffraction  grating  can  also  be  considered  as  an  infinite 
series  of  points  with  spacing  2a,  which  is  observed  through  a 
window  that  is  transparent  inside  the  limits  ±  w  and  opaque  out¬ 
side.  See  Figs.  17,  18.  Since  f(x)  is  the  product  of  two  functions 
the  diffraction  pattern  is  the  convolution  of  their  patterns 
according  to  (48) 

g  (u)  ^fgi  («)  ga  (u-d)  dd  (80) 

Now  the  pattern  for  the  infinite  grating  degenerates  to  the  series 
of  equally  spaced  points  shown  in  Pig.  18,  so  that  it  acts  as  a 
window  of  vertical  transparent  lines  of  spacing  u  n/a.  Thus,  to 
a  first  approximation,  the  pattern  is  the  single  slit  pattern 
traced  by  the  line  nearest  the  peak.  A.  L.  Patterson  has  applied 
this  point  of  view  to  the  diffraction  of  x-rays  by  small  crys¬ 
tals.  * 


J _ L! _ 1 

i  1 

1 

II  1  1  n  '  >  t  1  1 

-•H - h-2w  — M - h 

Fig.  17.  Orating  as  Product  of  Infinite 
Orating  and  Window  of  Width  2w 
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Fig.  18.  Orating  Pattern  aa  the  Convolution 
of  (eln  au)Ai  and  the  eerie*  of 
equidistant  lines  of  separation 
Av  •  au  ■  Tt 


6.  Grating  with  Finite  Slits  or  Dipoles 

Since  an  actual  diffraction  grating  consisting  of  n  slits  of  eidth 
2b  is  the  convolution  of  the  point  grating  f|^(x)  and  the  silt  of 
width  2b j  Its  diffraction  pattern  is  simply  the  product  of  (88)  by 
the  single  silt  pattern  (sin  bu/u) • 

8 («)  *  gn  bu  ^  b  sin  nau  sin  bu 

u  sin  au  bu 

A.t  the  peak  this  reduces  to  nb  which  Is  the  total  opening  of  the 
grating.  Obviously^  any  other  pattern  such  as  that  of  a  dipole  can 
be  substituted  for  the  pattern  from  a  single  opening.  In  any  ease 
It  cannot  effect  the  close  side  lobes  but  does  help  diminish  the 
higher  orders. 

0.  CONVOLUTIONS  OP  LXNB  80URCB8  IN  ONB  DIMBNSION 


In  the  following  exavples  the  convolution  usually  overlaps.  Pig. 
19  shows  how  a  uniform  line  source  of  width  8a  can  be  split  into 
a  number  of  seetlons«  each  of  which  is  spread  out  into  a  second 
uniform  line  source  of  width  Sb.  We  will  not  worry  much  about  the 
normalisation  constant  but  will  usually  1st  g(0)  ■  1  which  Implies 
that  the  area  of  f(x)  is  normalised  to  unity. 

1.  TrapeMoidal  Illuuination  (Pig.  19) 

The  result  of  the  ccmvolutlon  of  Pig.  19  is  a  trapesoldal  ampli¬ 
tude  distribution,  the  two  bases  being  8(b’^a)  and  8(b-a).  The 
diffraction  pattern  is  proportional  to  the  product  of  the  patterns 
of  the  two  line  sources. 


g(u)  e  gtgt 


sin  au  sin  bu 
au  bu 


(81) 


Since  the  numerator  cannot  exceed  unity,  the  amplitude  of  the  side 
lobes  cannot  exceed  a  fraction  l/(abu  )  of  the  peak. 


X*  1 
ab  (ansln  9)* 


(68) 
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Types  of  f(x) 


«(«) 


«i  (u)  - 
(«)  = 


Si  St 


sin  an  sin  bu 
an  Bu 


Figure  21 •  Gable 
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Figure  22«  Step  Approximation 

f„(x)  -  (1+1)*-  1,4,G,4,1 


a.  lu)  sin  an 

'  — Sii 

cos*au  sin  au 
au 


FIOHRBS  19  -  i69.  CONVOLUTIONS  OF  UNIFORM  LINE  SOURCES 


7BL4 


Thus,  the  intensity  decreases  as  (sin  0)  due  largely  to  the 
taper  at  the  ends*  Convolutions  of  three  line  SQurces  >111  de¬ 
crease  at  an  even  greater  rate. 

2.  Triangular  Illumination  (Fig.  20) 


When  a  -  b,  the  trapezoid  degenerates  to  the  triangle  of  Fig.  20 
with  an  overall  width  of  4a  and  a  pattern 


«(u) 


[ 


sin  au 
au 


■] 


2 


(63) 


3.  Gable  Illumination  (Fig.  21) 


The  gable  amplitude  illumination  Pig. ,21  is  the  sum  of  a  triangle 
of  hel^t  A  and  a  rectangle  of  hel^^t  B.  The  pattern  Is: 

A  teln  au J  *  +  B  sin  2  au 

2  [  au  y  2au  (64) 

4.  Approximation  of  t(x)*by  its  Value  at  n  Equidistant  Points 

If  an  amplitude  function  f(x)  Is  listed  at  n  equidistant  points 
^Be  diffraction  pattern  gf)(u)  corresponding  to  ^^(x)  elll 
show  strong  minima  and  maxima  analogous  to  the  orders  of  a  dif¬ 
fraction  grating.  In  order  to  keep  the  errors  small  the  phase 
between  adjacent  points  must  be  kept  small.  These  errors  are 
partially  corrected  in  the  methods  described  below. 

5.  Approximation  of  f(x)  by  Steps 

In  Fig.  22  the  anqilltude  function  f  (x)  Is  divided  into  rectangular 
steps  of  equal  width  2a  and  heights  f,,  equal  to  the  value  of  the 
ordinate  at  the  midpoints  of  the  Intervals.  The  resulting  approx¬ 
imation  Is  the  convolution  of  the  point  function  end  the 

uniform  line  source  of  width  2a,  so  that  the  approximate  pattern 
Is  the  product  of  the  transform  of  single  slit: 

gn(u)  sin  au 

au  (66). 


The  g||(u)  is  calculated  exactly  as  In  the  preceding  paragraph.  The 
added  factor  (sin  au)/au  Improves  the  approximation  by  suppressing 
the  orders  Introduced  by  the  grating,  f^(x).  Fig.  22  Is  drawn  for 
the  case  In  which  proportional  to  1,  4,  6,  4,  1,  the 
binomial  expansion  (l'*’l)  . 


6.  Approximation  of  f(x)  by  Chords 


Fig.  23  shows  how  a  function,  which  Is  zero  at  Its  end  points,  can 
be  approximated  by  a  series  of  connecting  chords.  These  are 
eqiuvalent  to  the  convolution  of  f|f(x)  and  a  triangle  of  width  4a. 
The  diffraction  pattern  is  the  product 


«n  ^8ln  auj 


(66) 
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Should  f (x)  have  end  points  not  zerOf  say  fi  and  fjj,  then  a 
straight  line  function  (A+Bx)  connecting  fj  and  can  first  be 
subtracted  from  f(x).  The  pattern  of  (A+Bx)  is  then  calculated 
Independently  and  added  to  the  solution  (66) • 

H.  BiroCT  OF  A  GAP  IN  A  LINS  80URCB 

Occasionally  an  array  Is  constructed  as  in  Fig.  M  with  a  narrow 
gap  at  the  center.  The  effect  of  the  gap  is  to  Increase  the  in¬ 
tensity  of  the  first  side  lobe  from  a  fraction  p*  to  a  fraction 
approximating 

(p+q)*  /(1-q)*  (67) 

where  q  is  the  fraction  of  the  area  of  the  amplitude  function  f(x) 
cut  out  by  the  gap. 


1^^ 

0 

•fO 

<1 

fgU) 

•mm 

Fig.  M  Gap  in  Array 

Ratio  of  Areas  is  q 


Fig.  M  Gap  Increases  First  Side 
Lobe  by  Raising  Axis 


This  la  due  to  the  fact  that  the  actual  amplitude  f(x)  of  the 
array  la  the  difference  of  the  original  function  fi  (x)  and  the 
gap  fc  (x) .  The  pattern  gt  (u)  for  the  gap  la  so  broad  that  it  has 
the  effect  of  elevating  the  axis  of  Fig.  SB  by  the  fraction  q. 
Since  the  amplitude  of  the  first  aide  lobe  was  originally  -p«  the 
new  side  lobe  has  a  relative  emolltude  - (p+q)/(l-q) ,  and  a  rel¬ 
ative  Intensity  of  (p*q)  /(1-q)  .  The  second  aide  lobe,  on  the 
contrary,  is  reduced  and  may  even  vanish. 


This  formula  holds  exactly  when  the  central  point  source  Is 
blocked  out  from  an  odd-numbered  n-polnt  diffraction  grating.  In 
which  case  q  =  1/n.  The  Intensity  of  the  side  lobes  midway  between 
two  orders  was  previously  1/n  .  The  odd  ones  are  now  Increased  to 
4/(n-l)  while  the  even  ones  vanish.  When  applied  to  a  3  point 
array  the  relative  intensity  of  the  one  and  only  side  lobe  In¬ 
creases  from  1/9  to  1. 

I.  INTEQRilTlON  BY  PARTS 


This  method  can  be  employed  when  there  are  discontinuities  in 
f  (x)  or  Its  derivatives  at  only  a  few  points,  Xj^.  Between  two 
such  points.  Integration  by  parts  yields: 

*»  -  ** 


,  ^  lux  lux 

f(x)e  dx  -  e 


{jL.  ^  ^  *  -DL_...7f(*) 

I  lu  (iu)»  (lu)*  I 

^  j  Xi 


On  summing  over  all  Intervals  the  contributions  cancel  except  for 
the  negative  difference  of  D'^f  (x)  ,  (x) ,  at  each  discontinuity. 

. (») 

We  have  so  far  Ignored  the  pattern  X^k*^'**  due  to  a  set  of 
radiant  points. 


F(x)  = 


f  (x)dx. 


then  each  radiant  point  can  be  considered  as  a  discontinuity  In 
F  (x) .  Since 


f(x)  .  DF(x) 


the  complete  formula  is 


where  xj^  is  now  the  position  of  either  a  radiant  point  or  a  dis¬ 
continuity  in  f(x).  Thus,  If  D^(x)  is  the  lowest  order  derivative 
containing  a  discontinuity,  the  leading  term  in  g(u)  will  vary  as 
l/u**  and  the  envelope  of  the  power  pattern  will  vary  as  1/u*". 
The  corresponding  condition  on  f(x)  is  that  D^f (x)  is  the  lowest 
order  derivative  which  becomes  infinite.  Usually  the  illumination 
and  Its  derivatives  are  continuous  except  at  the  ends  of  the 
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interval.  If  f(x)  Is  also  an  even  function  with  a  range  ±  1, 


g(<p) 


f  (1) 


+  cos  <1> 
cp 


(X> 

(72) 


Several  examples  of  Eqs.  (71)  and  (72)  are  given  in  Fig.  26 


f  (x) 


-0  +a 


g(u) 


lua  -iau 

e  +e  =2  cos  au 


Fig.  26  Examples  of  Diffraction  Patterns  Derived  by  Eqs.  (71,72) 

Examples  (1),  (2),  (3)  show  discontinuities  in  F(x)  and  Its  first 
and  second  derivatives  respectively.  (See  Eq.  71).  The  results 
In  (2)  and  (3)  could  also  be  obtained  from  (72) .  Examplesd  and  6 
combine  two  or  more  orders  of  discontinuities. 

J.  SraiBS  FOR  DIFFRACTION  PATTERN**'** 

For  moderately  small  values  of  u  it  is  feasible  to  expand  the  e^^* 
in  the  Fourier  integral  in  a  power  series  and  then  to  integrate 
term  by  term  as  follows. 

1.  Diffraction  Pattern  g(u)  in  Terms  of  Moments  of  f(x) 


26 
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g(u)  «C  Jr(x)  e^dx  .J*j^+lux  +  (lux) *  +...J  f  (x)dx 

“  Jf(x)dx  +  lurxf(x)dx  +...+  l"u^  |x"f(x)dx  ♦••• 

nl  J 

Let  us  designate  by  the  kth  oonent  of  f(x), 

Ujj  >  f  (x)  dx 


Then  (73)  becomes 


g(u)  «C  m  jit  ^  U4  u’ 

^  kl  '  31 

*  1  [>‘.<--■57  »•  ♦ 


An  even  function  has  only  even  moments  vhlle  an  odd  function  has 
only  odd  moments*  The  idea  of  a  moment  is  common  in  mechanics*  The 
zero  moment  is  the  area  /f(x)dx*  The  center  of  gravity,  c*g*, 
and  the  square  of  the  radius  of  gyration,  r‘,  of  a  distribution  of 
mass  points  are  also  simply  expressed  in  terms  of  momenta* 


/xf(x)dx  Ui 
/f(x)dx  u- 


/x*f(x)dx 

/f(x)dx 


The  series  (76)  converges  when  f(x)  Is  positive  and  confined  to 
the  range  ±  a,  for  then,  and  the  terms  of  the  series  do 
not  exceed  those  of  the  series  Poe^*”,  which  converges*  For  the 
case  when  f(x)  consists  of  two  points  of  amplitude  i  placed  at  ± 
a,  ■  a^  and  the  series  (75)  reduces  to  cos  au*  Again,  for 
uniform  illumination  between  *  a 


g(u)  ”  2a  )  1  >  -  ***  I  -  2a - 

L  81  51  J  V 

a.  Series  for  the  Power  Pattern  P(u)  ■  g(u)  g(u)* 

The  series  for  the  power  diffraction  pattern  is  the  sum  of  the 
squares  of  the  real  and  imaginary  series  of  (75)  according  to  (29) 

X  -  ..  «  / . «v.t  ^  Ug‘]  ..4 


>*o*  "  *f- 

M-  J.  ^  .  ' 


*.  /^O  —  /*»/>  V  ti  « 

■<■)  ■••••  ctmrln  eoajutat*  of  c(a) 


)  u 
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Figure  ZT •  Comparison  of  Approximations  Over  the  Main  Lobe  of: 

(a)  Two  Points,  (b)  Pnlfor*',  Illumination,  (c)  Half  Cycle  Cosine 
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Let  us  assume  that  Uj  is  zero,  which  is  true  for  an  even  function, 
or  for  any  function  whose  origin  is  shifted  to  its  center  of 
gravity,  (76).  On  factoring  out 


We  may  employ  the  series 

s 

N  =  -loggP(u)  =  -logg(l-z)  =  z  + —  +  ••• 


to  obtain  the  power  drop  N  in  Napiers. 

4  u,  2  r 

U  m-^  U  M 

“■>  I 


3‘ 


The  term  in  parenthesis  is  a  small  quantity  which  approaches  zero 
as  P(u)  approaches  the  Gauss  Error  curve  .  Th»  Hi>nn  n 


(81) 


The  drop  n 


In  db's  IS  4.343  times  N. 

Vtdths  of  Diffraction  Patterns 

The  series  method  can  be  employed  to  solve  for  the 
diffraction  patterns.  On  reversing  the  series  (81) 


p/z  D 

u\  E  -  k  4? 

widths  ^f  ^ 


u=  =  iio  nTi  .A.  JVL*)  1] 

U2  L  (  31^2  7 


iio  JT 

A  -  / 

-1 

If 

Us  * 

L  \ 

3U,V 

8  J 

where  A  and  B  are  constants.  Eq.  82  is  an  improvement  over  that  in 
RL  Report  T-7*  which  was  developed  as  follows.  If  P  is  the  frac¬ 
tional  drop  in  amplitude  measured  from  the  peak,  then  g(u)  In  (75) 
equals  (1-3) .  If  we  neglect  terms  above  and  solve  for  u,  then 


u 


(83) 


If  we  include  the  term  in  U4  but  assume  that  It  is  the  same  as  for 
the  Gauss  Error  curve,  then  the  following  approximation  results 
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(84) 


This  18  equivalent  to  Eqs.  (IS)  and  (98)  of  Report  T-7.  In  the 
following  discussion  we  use  cp  =  au« 

The  several  approximate  formulas  are  compared  in  Fig.  27  for  the 
main  lobe  of  the  patterns  for  three  illumination  functions: 
(a)  two  points,  (b)  uniform  illumination  and  (c)  the  half  cycle 
cosine.  The  simplest  approximation  is  cp^,  based  on  the  Oauss  Error 
curve.  This  is  fairly  good  down  to  half  power.  Eq.  (84),  indicated 
<P3*  lies  inside  by  a  fraction  which  depends  only  on  n. 
Although  q>^  may  fit  some  patterns  it  cannot  compare  with  the 
surprisingly  accurate  approximation  based  on  (82) .  The  percent 
error  in  width  using  is  shown  in  Fig.  28.  For  the  case  of 
uniform  illumination  the  error  is  less  than  0»b%  down  to  26  db. 
Paradoxically,  the  reversed  series  (82)  is  more  accurate  than  the 
direct  series  (81) ,  which  in  turn  is  much  closer  than  the  original 
arithmetic  senes  (79). 


-l%i 


20  db 


30  db 


Fig.  28.  Percent  Error  in  Width  Formula^  W  ==  A  iTnd-Bn) 

For  Three  Types  of  Illumination 

4.  Calculation  of  Beam  Widths 

The  been  width  W  in  radians  of  a  diffraction  pattern  is  approx> 
lasted  by 


W  ■  2  sincp 

whioh  IS  (X/r[}u.  Thus  we  may  write 


(86) 


30 


762-1 


w  -  A  fJTd-Bn) 


(86) 


where  X  and  the  diameter,  d  * 
constants  A  and  B  are  defined  as 

A  »  0.806  a\pi£  ,  B  - 
»  U« 


8a,  are  In  like  units,  and  the 
follows 


Table  1  Constants  For  Width  Formula 


Type  of 
f  (x) 

D  D 

Cos  h  nx 

A 

.306 

.629 

\ 

.702 

B 

.0198 

.012 

.0079 

n 

6 

26 

11 

Table  1  lists  the  constants  A  and  B  tor  use  in  (86)*  The  n  In  the 
table  is  the  range  in  dbs  over  which  the  error  In  width  is  less 
than  0.6)^. 


Bach  constant  is  dlmenaionlese.  For  example,  for  uniform  illumina¬ 
tion,  *  SaV  (k+l)  end 


W 


0.«89  fn 


(■ 


(88) 


It  is  an  interesting  coincidence  that  when  W**  is  In  degrees,  X  In 
centimeters  and  d  in  feet,  that  A  ■  0.996  or  approximately  unity 
Thus 


(uni form  llltBlnatlon) 


(89) 


5.  Beam  Widths  of  Convolutions 


The  pattern  of  a  convolution  is  the  product  of  too  patterns 
gi  St«  which  usually  approximates  a  Gauss  Error  curve  even  closer 
than  does  either  pattern  separately.  We  can  thus  drop  the  B  in 
(86) .  Since  the  second  moments  are  additive,  it  follows  that  the 
width  W  of  the  convolution  is  related  to  the  widths  W,  and  W,  of 
gt  and  gg  by  the  expression 


(90) 


768-1 


81 


A  Simple  graphical  solution  is  given  in  Fig 


29 


Fig.  29.  Relation  Between  Widths  W,  and  W^. 


The  author  wishes  to  thank  Uary  Batchelder  for  the  calculations  in 
Figs.  27  and  28,  and  Pauline  Austin  for  her  careful  reading  of 
the  manuscript. 


R.  C,  Spencer 
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